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Abstract. We consider a spin-boson model in which a spin 1 system is coupled to 
an oscillator. A unitary transformation is applied which allows a separation of terms 
responsible for the Bloch-Siegert shift, and terms responsible for the level splittings at 
anticrossings associated with Bloch-Siegert resonances. When the oscillator is highly 
excited, the system can maintain resonance for sequential multiphoton transitions. At 
lower levels of excitation, resonance cannot be maintained because energy exchange 
with the oscillator changes the level shift. An estimate for the critical excitation level 
of the oscillator is developed. 
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1. Introduction 

The interaction of a two-level system with oscillatory off-diagonal coupling leads to a 
shift in the transition energy, known in the literature as the Bloch-Siegert shift [TJ [2]. 
When the shifted transition energy is resonant with an odd number of oscillator quanta, 
energy exchange between the two systems can occur. This effect appears in energy level 
calculations as a level splitting at the Bloch-Siegert resonances. Both the shift and 
splitting have been studied previously using models based on the Rabi Hamiltonian (in 
which the oscillator is presumed classical) [H [21 [31 HI El [6] and based on the spin-boson 
Hamiltonian (in which the oscillator is presumed quantum mechanical) IB]- It has 
been noted that the Bloch-Siegert shift arises from the counter-rotating terms in the 
Hamiltonian, and that it disappears when the rotating-wave approximation is made 

mu ■ 

Our interest in this problem is motivated by the possibility of coherent energy 
exchange between quantum systems with highly mismatched characteristic energies. For 
example, we have been interested in the dynamics of energy exchange between a two- 
level system with a large transition energy AE, and an oscillator with a small energy 
quantum Hujq. The spin-boson Hamiltonian is one of the simplest models exhibiting 
such coherent energy exchange. We recently considered the Bloch-Siegert shift [10] and 
level splittings [H] at the Bloch-Siegert resonances, the latter of which is a result of the 
coherent energy exchange between the two-level system and oscillator. 

Most of the work cited above is focussed on the two-level problem, which is usually 
modeled as a spin 1/2 system interacting with a single harmonic oscillator. The 
generalization of the problem to systems involving higher spin leads to more complicated 
models which have not received comparable attention in the literature. The Bloch- 
Siegert shift for the spin 1 case was studied previously by Hermann and Swain [12], [T3] . 
Some progress has also been made in the case of the general spin problem [B]. 

In this work we focus on the problem of a spin 1 system coupled to a simple 
harmonic oscillator. Our earlier analysis in the spin 1/2 version of the problem made 
use of a unitary transform in which the shift and level splitting effects can be associated 
separately with different terms in the rotated Hamiltonian. Such a separation serves as 
the basis for the development of analytic results for both shift and splitting which are 
useful over a wider range of coupling strength than available previously. The unitary 
transform that we used for the spin 1/2 problem can be extended simply to higher spin 
models, allowing for a separation of shift and level splitting effects in more complicated 



Multiphoton Bloch-Siegert shifts and level-splittings in spin-one systems 



3 



problems. We have decided to focus here on the spin 1 case since it provides a good 
example of this generalization. 

Although the spin 1 system is a three-state system, we have found it useful to think 
about it in terms of an underlying two-spin problem. For example, the Bloch-Siegert 
shift that we find below for the spin 1 can be understood simply as arising from the 
individual shifts associated with two spin 1/2 systems. Later on in this work, we find 
that the Bloch-Siegert resonance condition cannot be maintained at modest n, which 
can be understood in terms of an initial resonant spin 1/2 transition that exchanges 
energy with the oscillator, followed by a second spin 1/2 transition that is no longer 
resonant since the oscillator is changed. 

2. Unitary equivalent Hamiltonian 

We focus on a spin 1 system coupled to an oscillator, using a spin-boson Hamiltonian 
of the form 



We are interested in the regime where the photon excitation n is large, and where the 
transition energy is much greater than an oscillator quantum (AE 3> Huq). 

As was the case in the spin 1/2 problem [TUl [TTj . it is useful to consider in the case 
of the spin 1 problem the unitary equivalent Hamiltonian 




(1) 



H' 



(2) 



where 




The rotated Hamiltonian becomes 



H' 



H + V + W 



(3) 



where 




(4) 
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Within the parameter space of interest to us, the term is small; consequently, we 
neglect it in what follows. 



3. Approximate energy eigenvalues 

In the large n limit, the oscillator has a strong impact on the spin system, but the spin 
system impacts the oscillator only weakly. As a result, the energy levels are found to 
be reasonably well approximated away from the level anticrossings by 



E n>M {g) = AE(g)M + huj n 



(7) 



where AE(g) is the dressed two-level system energy. This is very much like the behaviour 
found for the spin 1/2 version of the problem discussed previously [10J. It is useful to 
adopt the same definition for the dimensionless coupling constant g for the spin 1 version 
of the problem 

UJn^ 



AE 

that we used before for the spin 1/2 case. 

The dressed two-level system energy AE(g) is also the same, as would be expected 
since it comes about from the same basic interaction. We can think of the spin 1 system 
as being made up of two spin 1/2 systems coupled identically to the same oscillator. 
Since the Bloch-Siegert shift is the result of the interaction between a single two-level 
system and oscillator, we should expect that it will be nearly the same per two-level 
system if there are more than one. 



Multiphoton Bloch-Siegert shifts and level-splittings in spin- one systems 5 

We can make use of the rotation to see this. Consider eigenf unctions of the 
unperturbed Hq part of the rotated Hamiltonian 

E nt M?Pn,M = Holft n ,M (8) 

We can separate variables (as done in the spin 1/2 case) to write 

V>n,M = \S,M)u n>M (y) 
where m„,m(2/) satisfies 

u n ,M + M^AE 2 + 8U 2 y 2 u nM (9) 

In the large n limit, this is just a modified simple harmonic oscillator with a small in- 
dependent perturbation. By adopting a simple harmonic oscillator wavefunction as a 
trial solution in a variational computation, we obtain approximate energy levels of the 
form of Equation (J7J) where 

AE(g) = L 

This dressed transition energy, identical to what we found in the spin 1/2 case [TU] , is 
in good agreement with the numerical results away from the level anticrossings. 

4. Multiphoton resonances 

As mentioned above, our interest in this model is driven in part by the possibility 
of coherent energy exchange between systems with strongly mismatched characteristic 
energies. In the spin 1 problem, resonances occur associated with anticrossings when 
the dressed two-level system energy AE(g) matches an odd number of oscillator quanta. 
In the large n limit, the oscillator is impacted only weakly through the exchange of 
a modest number of oscillator quanta; in which case the dressed energy AE(g) is 
approximately invariant, and we are able to develop a resonance condition applicable 
to both transitions. When n is smaller, the change in the number of oscillator quanta 
may be a significant fraction of the total number of oscillator quanta; in which case 
the dressed two-level transition energy may be on resonance for one transition and off 
resonance for the other. We begin our discussion here focusing on the first case, since 
it is simpler. 



E n ,M + 



U n ,M 



HuJn 
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(10) 
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4-1. Large n resonance condition 

We are interested then in the resonance condition at which the dressed transition energy 
is matched to an odd number of oscillator quanta. To proceed, we consider three basis 
states which are eigenstates of the rotated Hamiltonian H : 

<j)-i = u n+ An,-i(y)\S,-l) (po = u nfi (y)\S,0) 0i = Un-A. n) i(y)\S, 1)(11) 

where S = 1. These states have energies e_i, eo and e% respectively, which in the large 
n limit are given approximately by 

e_i = huj (n + An) - AE{g) (12) 

e = hu Q n (13) 

ei = hu (n - An) + AE(g) (14) 

with AE{g) taken to be a constant for both transitions. The resonance condition in 
this case is 

AE(g) = An?iuj (15) 

where, as in the two-level case, the number of oscillator quanta exchanged An must be 
odd for level splitting to occur. 

4-2. Comparison with previous work 

As mentioned above, Hermann and Swain [13] have calculated the Bloch-Siegert shift 
for the spin 1 case. In their notation, the resonance condition to fourth order is 

m ujq 2V 2 p p 3 (3p 2 -7)F 4 

p (p 2 - iH (p 2 - i) 3 ^o 3 1 ' 

In order to write this resonance condition in our notation, we make the following 
replacements 

V -> UV2n~, u min -> uo, uj — > AE, p^ {2k + 1) 

to obtain 



AE 



i 4 g 2 (2fc + l) 2 _ 4(2fc + l) 4 [3(2fc + l) 2 - 7] g 4 
+ 4k(k + l) ~ [4A;(fc + l)] 3 + 



(2k+l)hu (17) 



Note that Equation ( ITTj) is exactly the same condition as obtained by Ahmad and 
Bullough [21 [10] for the spin 1/2 problem. 
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To compare our calculation to these perturbative results, we can expand AE(g) 
from Equation ( JTU1) in powers of g to obtain the condition for the (2k + l)th resonance 
as 



where An = 2k + l. This condition is the same as was found in the spin 1/2 problem [10J. 
In the limit that k 3> 1, we see that Equation (II7j) reduces to our result [Equation (fT8l) ]. 

5. Dynamics: resonant case 

The basis states that we chose [Equation (TTT1) ] are matched total energy states; the 
spin 1 system pieces are combined with modified oscillator states in which an odd 
number of oscillator quanta are matched to the dressed two-level transition energy. 
When the coupled system makes transitions between these different basis states, the 
resulting dynamics describe coherent energy exchange between component quantum 
systems with very different characteristic energies. The resonant case in the high n 
limit is most interesting in this regard, since both two-level systems make transitions, 
and a full 2An oscillator quanta are exchanged. 

5.1. Three-state model 

To proceed, we consider a dynamical state ip(t) constructed from the three basis states 
listed in Equation (TTTT) 



In the large n limit, these three basis states become degenerate. In addition, the different 
matrix elements become roughly equal 





C_i(*)0_i + C O (£)0o + Ci(t)(f) 



1 



(19) 




(20) 



where 



v = 



2huj U 
AE 



(21) 



and 




(22) 
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Figure 1. Occupation probabilities as a function of time for the degenerate three 
state model. The time axis is in units of uj$t where ujq = V2v/h. 



' d_ 

dy 



U n +An,Al(y) 



1 



\AEj 



>u n ,M(y)dy 



(23) 



The two integrals appearing here have been found to approach a common limit (to 
within the sign) based on calculations using numerical solutions for the rotated frame 
Hq problem, and based on calculations using the WKB approximation. This behaviour 
can be understood simply by noting that the the derivative can be expressed in terms 
of raising and lowering operators and in the large n limit the u n (y) functions are very 
nearly pure harmonic oscillator states. 

The inclusion of the perturbation V leads to the restricted Schrodinger equation 

/ c_i \ / e v \ / c_j \ 



in— 
dt 



Co 

V Ci J 



v e v 
V v e J 



Co 

V c i J 



(24) 



5.2. Dynamical solution 

Solutions for Equation (1241) can be readily constructed using an eigenfunction expansion. 
One interesting solution is 



c-i(t) 



co(t) 



( 1 

2 



-ieot/h 



1 

V 2 



COS + 1 



_j_ ■ V2vt 



cos ^ - 1 



(25) 
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Figure 2. Expectation value (M) and (n — jiq) as a function of time. The time 
axis is in units of uj^t where ujq — V^v/h. This computation is for a resonant energy 
exchange process where 25 oscillator quanta are exchanged for one dressed two-level 
system quantum. 

The associated probabilities are illustrated in Figure [TJ One sees that the system starts 
in state 0_i with unity probability, moves through state </>q, and then reaches <\>\ with 
probability unity. The solutions are periodic, so that the system cycles back and forth 
between the different states. 

5.3. Evolution of expectation values 

In Figure [2] we show results for the expectation values (M) and (n — no) in a resonant 
energy exchange process where 25 oscillator quanta are exchanged for a single dressed 
two-level system quantum. The expectation values in this case are computed according 
to 

{M)=Y J M'\c M '(t)\ 2 (n-n ) = - £ M' An\c M > {t)\ 2 (26) 

M' M> 

One sees in this case a complete energy exchange between the two systems, where the 
excitation energy of two dressed two-level systems [2AE(g)\ is exchanged for an equal 
amount of oscillator energy (2AntiLdo). 

5.4- Oscillations of (M) are sinusoidal 

We can look at the dynamics of energy exchange in the large n resonance case in 
another way. Within the set of basis states that we have selected, the basis states 
are degenerate on resonance, and the coupling is proportional to matrix elements of S x 
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[see Equation (Ell) ]. In this case, we can determine the evolution of (S z ) for the spin 1 
case (or in other cases as well) through commutation with the interaction term of the 
associated restricted Hamiltonian. The result is 

^(S z ) = -tf(S z ) (27) 

where the frequency Q is 

Hence if all the two-level systems are initially in the ground state (in the rotated frame), 
in general the average (S z ) will exhibit sinusoidal dynamics according to 

(S z ) = (M) = -Scos{Qt) (29) 

This is consistent with the results shown in Figure 2. Note that these dynamics are in the 
rotated frame, and could be obtained if the interaction U were turned on adiabatically. 

6. Finite n effects 

Up to this point we have assumed that n is "large enough" so that the basis states 
are degenerate. However, when n is not overly large, the levels are not degenerate and 
hence it is not possible to arrange for a clean energy exchange as described above. This 
motivates an interest in understanding how large n must be so that the system acts as 
if it is degenerate. 

6.1. Parameterization of basis state energies 

For finite n, we have found (based on numerical calculations) that the energy levels of the 
original Hamiltonian H can be accurately fit away from resonance using the expansion 
of the form 

E nM = A + B(n - no) + CM + DM(n - n ) + FM 2 + ■■■ (30) 

where A, B, C, D, and F are fitting coefficients. The terms quadratic in An are not 
that significant because the oscillator is highly excited. We can relate the basis state 
energies e_i and e\ to eo using this parameterization to obtain 

e _! = eo + BAn-C - DAn + F + --- (31) 

€l = eo -BAn + C-DAn + F + --- (32) 
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We take the approximate resonance condition to be 
BAn = C 

which is nearly equivalent to our resonance condition from above 
Antuvo = AE(g) 
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(33) 
(34) 



6.2. Mismatch in basis state energies at resonance 

When this resonance condition is satisfied, the energy levels (in the absence of coupling 
terms) are still not matched. Instead, we find that e lies above e\ and e_i due to the 
presence of higher-order terms in the fitting expansion 

DAn + F + ■ ■ ■ (35) 



£-1 — eo 
6l = e - DAn + F H 

Both analytic and numeric computations lead to the conclusion that 

,_, ^ . constant 
F - DAn 



n 



(36) 



(37) 



for large no- The numerator on the RHS is a constant that can be determined from 
parameterizing the levels in a direct numerical calculation, or can be estimated from 
perturbation theory in the rotated frame. For the purposes of discussion, we may write 
this as 

[no(F-DAn)} 

no=oo 



F- DAn 



n 



(38) 



6.3. Determination of critical n where splitting matches coupling 



The matrix element that produces transitions in the rotated frame can be approximated 
by 



{<Po\V\^ 



2huoU 
~AE~' 



2hu g 



(39) 



UQ=CO 



We can now estimate the number of oscillator quanta n required to make the coupling 

matrix element equal in magnitude to the basis state splitting by requiring 

[n (F - DAn)] n0=oo 



\E\ — Eq\ 



n 



y/2\((t>o\V\(t>i) 



(40) 



The y/2 that appears here reflects the extra factor in the level splittings obtained from a 
diagonalization of the three-state model on resonance. This is satisfied when n is equal 
to a critical value 

[n (F - DAn)] n0=oo 



Tl>crit 



2V2hu} g I/y/n~o~ 



no=co 



(41) 
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Figure 3. Level splitting from direct solution of original Hamiltonian H problem as 
a function of no. The level splitting is between states that in the rotated frame are 
mostly composed of M — — 1, An — 15, and M = 1, An = — 15, for a model in which 
AE = llSwo. 




6.4- Numerical example 

In Figure 3 is shown the level splittings for state computed from the original Hamiltonian 
H without rotation in an example that illustrates this effect. In this calculation, we have 
selected a model with AE = 11 hwo, and we focus on the resonance at AE(g) = 15 huiQ, 
where g is determined from numerical optimization to minimize the level splitting at each 
Hq. The numerical data shows a resonance behaviour similar to what we would expect 
from the three-state model above. At large no, the basis states are nearly degenerate 
relative to the coupling matrix element, and we find maximum level splitting. At small 
n , the basis states are separated by more than the coupling matrix element, so that 
the second-order coupling between 4>-i and 4>i is small, and hence the level splitting is 
also small. The critical number of oscillator quanta in this case is about 2.7 x 10 4 . 

6.5. Estimation of n crit from analytic estimate 

For this problem, the quantity I j \fn is found from a WKB calculation to be 

= 1.77 x 1(T 3 (42) 



n 

-l no=oo 

The D term estimated from 



D = (sSm e ^) m „ < 43 > 

\ / n=no,M=0 
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is used to obtain 

n D -> 6.60 QW) (44) 
as no goes to oo. From a direct parameterization of the energy levels we obtain 

n F -> -10.4 Qna; ) (45) 
We can compare this result from the value obtained using the WKB approximation 

n " F = T = - 1L0 (iW (46) 



These parameters are combined to produce the estimate based on Equation (jUJ 

n crit = 3.26 x 10 4 (47) 

This estimate is in reasonable agreement with the calculation of Figure 3. We conclude 
that a simple three-state model provides a good foundation for estimating how large n 
must be in order for the basis states to be degenerate relative to the coupling. 



7. Discussion 



We have previously applied a rotation in the case of the spin 1/2 problem in order 
to produce a dressed system in which the unperturbed Hamiltonian H provides a 
good approximation away from the level anticrossings at resonance [ID] , and where 
the perturbation V gives most of the coupling responsible for the level splittings at the 
anticrossings Here we have applied a similar rotation to the spin 1 problem, in which 
case a similar separation of the rotated Hamiltonian occurs. Although our discussion 
here is focused on the spin 1 problem, we have found that similar good results are also 
obtained in the case of higher spin as well. We find that the Bloch-Siegert shift in the 
spin 1 case is close to that of the spin 1/2 problem, in agreement with the perturbative 
result of Hermann and Swain [13] . 

We have considered level splitting at the anticrossing in association with the system 
dynamics at resonance. In the large n limit, energy exchange between the spin 1 system 
and oscillator constitutes only a small perturbation to the oscillator, so that the dressed 
transition energy is not changed, and that the same resonance condition applies to both 
transitions. In this case, it is possible for complete energy exchange to occur between the 
spin 1 system and the oscillator. We have given an analytic solution for the dynamics 
in the case of resonance. 
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At more modest n, the oscillator is modified sufficiently by a change of An oscillator 
quanta so that the resonance condition no longer holds for a second transition. In 
this case, we might think of the coupled system as being made up of two spin 1/2 
systems with different dressed transition energies weakly coupled to the oscillator. 
Accordingly, at larger n we might think of the coupled system as a dressed spin 1 
system weakly coupled to the oscillator. The transition between these two kinds of 
behaviour is determined by a critical excitation of the oscillator. We have developed an 
analytic estimate for this critical number of oscillator quanta n cr u [see Equation ( 14TT) ] 
in terms of fitting parameters which we can derive directly from solutions of the 
unrotated Hamiltonian H, from numerical solutions of the unperturbed part of the 
rotated Hamiltonian Hq, or from the WKB approximation as applied to the rotated 
problem. 

References 

[I] Bloch F and Siegert A 1940 Phys. Rev. 57 522 
[2] Shirley J 1965 Phys. Rev. 138, B979 

[3] Ahmad F and Bullough R K 1974 J. Phys. B: At. Mol. Phys. 7 L275 
[4] Fregenal D et al 2004 Phys. Rev. A 69 031401 (R) 
[5] F0rre M 2004 Phys. Rev. A 70 013406 

[6] Ostrovsky V N and Horsdal-Pedersen E 2004 Phys. Rev. A 70 033413 

[7] Cohen-Tannoudji C, Dupont-Roc J, and Fabre C 1973 J. Phys. B: At. Mol. Phys. 6 L214 

[8] Klimov A B, Sainz I, and Chumakov S M Phys. Rev. A 68 063811 

[9] Hattori T and Kobayashi T Phys. Rev. A 35 2733 

[10] Hagelstein P L and Chaudhary I U Preprint quant-ph/0709.1958 

[II] Hagelstein P L and Chaudhary I U Preprint quant-ph/0709.1961 
[12] Hermann J and Swain S 1977 J. Phys. B: At. Mol. Phys. 10 1815 
[13] Hermann J and Swain S 1977 J. Phys. B: At. Mol. Phys. 10 2111 



